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ABSTRACT 

We simulate the dynamics of the interaction between relativistically expanding spherically symmetric shells 
using a 1 -dimensional hydrodynamic code and calculate spectra and light-curves arising from such collisions by 
integrating the synchrotron and inverse Compton emission of the shocked gas. The numerical results reflect the 
most important features observed in Gamma-Ray Bursts: the spectrum exhibits a progressive softening (its break 
energy decaying exponentially with the 50-300 keV photon fluence), and the pulses that form the burst appear 
narrower in higher energy bands. Analytical results for the most important physical parameters of the burst are 
obtained by solving the shock jump conditions for a pair of interacting shells, in the case when both the forward 
and reverse shocks are relativistic. 

Subject headings: gamma-rays: bursts - methods: numerical - radiation mechanisms: non-thermal 



1. INTRODUCTION 

We consider Gamma-Ray Bursts (GRBs) arising from multi- 
ple internal collisions between relativistic shells of cold ejecta 
expanding into a vacuum. Under typical conditions, such col- 
lisions take place well before the decelerating effect caused by 
an external medium becomes important. The dynamics and en- 
ergetics of the unsteady outflow leading to internal shocks in an 
ejecta released, e.g. in a compact merger or collapse event, were 
outlined in Rees & Meszaros (1994). In this model, the physical 
conditions determining the energy deposition in the ejecta are 
not steady during the entire event, resulting in a non-uniform 
distribution of the bulk Lorentz factors within the ejecta. Faster 
shells of material catch up with slower ones, leading to the for- 
mation of two shocks in each pair of interacting shells. The 
slower shell is swept up by a forward shock that accelerates 
it, while the faster ejecta are decelerated by a reverse shock. 
Both shocks heat the ejecta, give rise to a turbulent magnetic 
field and accelerate electrons to a power-law distribution. The 
shocked fluid cools through synchrotron and inverse Compton 
emission and through adiabatic losses, due to the radial expan- 
sion of the ejecta. The radiation received by the observer comes 
from a small part of the ejecta, that moves almost toward the 
observer, due to the high Lorentz factor of the emitting fluid. 
Consequently, the burst seen by the observer is substantially 
blue-shifted and much shorter than as seen in the co-moving 
frame of the ejecta. 

The average spectra arising from internals shocks in such un- 
steady winds have been calculated analytically by Papathanas- 
siou & Meszaros (1996) for wide ranges of model parameters 
(magnetic field and electron acceleration efficiency, wind dura- 
tion and wind variability timescale, average Lorentz factor of 
the ejecta). Their results show that burst spectra can extend 
over many orders of magnitude, from eV to TeV. The spectra 
calculated by Pilla & Loeb (1998) show a similar wide range of 
photon energies, and a time evolution characterized by a soft- 
ening of the spectrum, determined by the electron synchrotron 
and inverse Compton cooling in the early phase of the burst, and 
dominated by the pair-cascade process at later times. For reason- 



able values of the outflow Lorentz factors and the shock radius, 
the depletion of high energy photons by pair creation was found 
to alter the burst spectra at energies higher than the BATSE win- 
dow. The two radiating processes which cool the electrons have 
been discussed by Sari & Piran (1997), who derived constraints 
on model parameters from the required efficiency, observed 
peak of the spectrum and ratio of the variability timescale to 
the burst duration and from the condition of optical thinness to 
pair creation. A detailed parameter search of the properties of 
time-integrated internal (as well as external) shocks provides 
constraints on the values of parameters which lead to bursts in 
the BATSE window (Papathanassiou & Meszaros 1998). Us- 
ing a code based on the kinematics of the interaction between 
many shells, Kobayashi, Piran and Sari (1997) calculated the 
efficiency with which shell collisions convert their bulk kinetic 
energy into internal energy and, assuming that this is radiated 
away, obtained a set of bolometric light-curves. In a more 
elaborate treatment, Daigne & Mochkovitch (1998) studied the 
temporal and spectral features of the synchrotron emission from 
a population of electrons altered by the power-law injection at 
shocks. They found a good agreement between most of their 
results and the features shown by real bursts: duration-hardness 
anti-correlation (e.g. Kouveliotouet al. 1993), appropriate tem- 
poral asymmetry of pulses (Norris et al. 1996), spectra well fit 
by the Band function (Band et al. 1993), spectral hardening be- 
fore a countrate increase (e.g. Bhatetal. 1995), pulse duration - 
energy dependence consistent with observations (e.g. Norris et 
al. 1996), and an exponential dependence of pulses peak energy 
on the photon fluence (Liang & Kargatis 1996). Other features 
of the real bursts, e.g. the general softening of the spectrum 
(Ford et al. 1995), are not well reproduced by the results pre- 
sented by Daigne & Mochkovitch (1998), probably due to the 
fact that they did not take into account the effects arising from 
the geometrical curvature of the emitting source or the radiative 
(synchrotron and inverse Compton) cooling of the electrons. 

In this work we investigate some of the most important prop- 
erties of internal shock bursts, using a 1 -dimensional hydro- 
dynamic code suitable for simulating the interaction between 
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perfect fluids, and calculate the light-curve and spectrum of the 
radiation in the BATSE range emitted by shock-accelerated rel- 
ativistic electrons, assuming spherical symmetry. In the frame- 
work of numerical hydrodynamic calculations we can follow 
in detail the post-shock time evolution of the electrons due to 
radiative and adiabatic losses in each cell of the fluid, and inte- 
grate the synchrotron output to determine the photon field that 
is up-scattered by the electrons. We then integrate the emission 
of the shocked fluid taking into account the dependence of the 
Doppler boosting on the angle between the radial direction of 
outflow and the direction toward the observer, and the effect 
of the shells' curvature on the photon arrival time. The bursts 
obtained with this time-dependent radiation and hydrodynamic 
treatment of the shock evolution show features that cannot be 
accounted for by kinematic treatments, including a softening 
of the spectrum and other correlations typical of the observed 
bursts. 



2. ANALYTICAL TREATMENT OF THE TWO-SHELL INTERACTION 

Before proceeding to a numerical hydrodynamical treatment, 
it is useful to start with some preliminary analytical insights. 
The physical conditions inside the shocked shells and the prop- 
erties of the emitted radiation can be estimated by first calculat- 
ing the Lorentz factor of the shock that sweeps up a shell during 
the interaction with another relativistic shell, in the frame of 
the yet un-shocked part of the shell. This can be done using 
the shock equations given by Blandford & McKee (1976) and 
the fact that the flow velocity and pressure in the shocked fluid 
has the same value on both sides of the contact discontinuity 
that separates the interacting shells. In the following analytic 
treatment, the internal pressure of the un-shocked fluid is ne- 
glected, therefore it applies only to a pair of cold shells. The 
effect of the radial expansion during the shells' collision is also 
neglected, therefore the following results are accurate only for 
shells so thin that they do not expand significantly before the 
shocks sweep them. If the two shells expand prior to their in- 
teraction in the way prescribed by Meszaros , Laguna & Rees 
(1993) (see below) then, during the time it takes the shocks to 
sweep up the ejecta, the radius of shells roughly doubles, there- 
fore the radial expansion cannot be neglected and the following 
analytic results should be regarded only as approximative. 

The inner shell has a co-moving frame density p, and moves 
at a Lorentz factor 1^, while the outer shell has a co-moving 
density p Q and moves at a Lorentz factor r o , lower than that 
of the inner shell. The collision of the two shells generates 
two shocks that compress and heat them. If Y s h denotes the 
Lorentz factor of the shocked fluid in the lab-frame, then the 
Lorentz factors Y' sh . and F' sh of this fluid in the frame of the 
yet un-shocked parts of the inner and, respectively, outer shell 
are: 
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The pressure of the shocked fluid around the contact disconti- 
nuity can be calculated for each shock: 



(psh)i/o = {Kh,i/o - l )^Kh,i/o + l )Pi/ c 2 



(4) 



where 7 is the adiabatic index of the heated fluid, assumed to 
be the same on both sides of the contact discontinuity. Using 
the equality of these pressures and equation (||), one obtains a 
quartic equation for the lab-frame Lorentz factor of the shocked 
fluid: 

7 (ff 4 - Y)x 4 + 2(7 - l)g(Y - .9 V +2(2 - j)g 2 (Y - l)x 2 + 



2tt-l)g(g 2 Y-l)x + *f(l-g 4 Y)=0 



(5) 



where Y = Pi/p a - 

Generally, equation ^) can be solved only numerically. 
Figure 1 shows the Lorentz factors r s h,i and r s h, of the 
shocked fluid as functions of the ratio of the pre-shock co- 
moving densities for three values of the ratio of the pre- 
shock lab frame Lorentz factors. For given parameter g, the 
T' sh i and T' sh curves are symmetric relative to the ordinate, 

i.e. T' sh 4 (Y) = F' sh (^ _1 ), as it can be shown using equations 
(||) and (0). It can be seen that, unless Ti/F <; 5 and either 
Y > lO^or Y- 1 > 10 2 , then T' sh < 2, which implies that 
the two shocks that sweep the shells (the reverse shock in the 
inner shell and the forward shock in the outer one) are only 
mildly relativistic. Figure 1 also shows the efficiency e at which 
the two shocks convert the shells' kinetic energy into internal, 
defined as the ratio between the lab-frame internal energy of the 
shocked gas and the kinetic energy that this gas had before it was 
swept up by one of the shocks. For given parameters g and Y, 
the efficiency is constant as long as there is a shock sweeping up 
each shell. The efficiency of kinetic to internal energy conver- 
sion changes after one of the shells has been entirely swept up 
by one of the shocks; in this case one can calculate the overall 
efficiency of the interaction from conservation of momentum 
and energy. The result is: 
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where p = Mi/M is the ratio of the shells' masses and 
G = g 2 = Ti /T is the ratio of their Lorentz factors. 

Analytic solutions of equation (Q) can be obtained in the 
particular case Y = 1 or if one makes the assumption that 
the two shocks are either quasi-newtonian or relativistic, which 
requires that Ti > r o orFi ^> r o , respectively. These solutions 
are: 
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where y — \/Y — \J pi/p - Once the shock Lorentz factors 
are known, one can calculate the co-moving internal and rest 
mass energy density of the shocked fluid, the turbulent mag- 
netic field B, and the minimum electron Lorentz factor j m of 
the power-law distribution in each shocked shell, from where 
all the important characteristics of the synchrotron and inverse 
Compton emission can be derived. Some analytical results in 
this direction are presented in Appendix A, only for the case of 
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relativistic shocks (T^ 3> r o ), as this case leads to higher burst 
efficiency and, thus, is more likely to be encountered in the 
bursts we observe. For definiteness we consider here the case 
of two shells that have not undergone any previous collisions, in 
which case the two unknown pre-shock co-moving densities pi 
and p can be correlated with basic burst parameters. The shells 
expand prior to the their interaction as described by Meszaios et 
al. (1993): they go through an acceleration phase, coast at con- 
stant Lorentz factor, and later start expanding, their lab-frame 
thicknesses evolving as r/I?, , where r is the radial coordinate. 
If the slower shell was released a time t v before the faster one, 
then the collision takes place at the "interaction radius" 

where the co-moving densities are 

ft/oC 2 = , (9) 

Ei and E Q being the kinetic energies of the two shells, therefore 

Y = E,/E . 

Equations (A6) and (A7) show that the burst spectrum de- 
pends strongly on the Lorentz factors of the two shells, the ratio 
of their energies, the burst variability timescale, and the electron 
acceleration efficiency. The dependence on the Lorentz factors 
and variability timescale arises mainly from the fact that we 
considered the interaction between shells that have propagated 
unperturbed until their collision; the results would be different 
if one or both shells interacted before with other shells, in which 
case equation ([)]) does not hold. The spectrum dependence on 
the ratios of the shells' Lorentz factors and of their energies 
comes from the hydrodynamics of the interaction. 



3. DESCRIPTION OF THE NUMERICAL CODE 

The numerical code that we developed contains two major 
parts: one which simulates the hydrodynamics of the interac- 
tion between two relativistically expanding fluids, described by 
Wen, Panaitescu & Laguna (1997), and one which calculates 
the emission of radiation from the shocked gases, through syn- 
chrotron and inverse Compton processes, and computes the ob- 
served spectrum and photon/energy light-curves, described by 
Panaitescu & Meszaros (1998a,b). The most important features 
of our calculation of the burst emission are listed below. 

Electron distribution and magnetic field intensity. The in- 
tensity B of the turbulent magnetic field is parameterized by the 
fraction e mag of the internal energy that is stored in the magnetic 
field. An electron power-law distribution (of exponent —p) is 
initialized in a grid cell containing shocked fluid when that cell is 
added to the shocked structure. The addition of a new "shocked" 
cell is done by the Glimm method that we use to simulate the 
propagation of the shocks, when the mass of the pre-shock fluid 
swept-up since the last added cell reaches the mass correspond- 
ing to the cell volume and post-shock density determined by the 
shock jump equations. The minimum electron random Lorentz 
factor 7 m of the power-law distribution is set by the energy given 
to leptons after shock heating, taken as a constant fraction e e i 
of the total internal energy of the newly shocked fluid, and by 
the fraction £ of electrons that are picked up by shock accelera- 
tion. The electrons are considered decoupled from protons and 



magnetic fields after their initial acceleration (i.e. they are not 
re-energized), and lose energy through emission of synchrotron 
and inverse Compton radiation and through adiabatic cooling. 

Radiative losses. Given the local value of the turbulent 
magnetic field B and the evolving electron distribution in each 
shocked cell, at some lab-frame time t, we calculate the syn- 
chrotron losses and integrate the emitted radiation over the entire 
volume of the shocked fluid and over the electron distribution, 
to calculate the synchrotron radiation energy density at each 
point in the shocked fluid, necessary for the computation of the 
inverse Compton losses. The approximations used for a faster 
numerical calculation of the synchrotron spectrum and inverse 
Compton losses are given in Appendix B. Also to reduce the 
computational effort, we do not use the full shape of the syn- 
chrotron spectrum when calculating the inverse Compton losses; 
instead we approximate as monochromatic the synchrotron ra- 
diation to be up-scattered, at an intensity-weighted frequency. 
Further, the spectrum of the up-scattered radiation is approxi- 
mated as monochromatic, at the peak frequency of the inverse 
Compton spectrum corresponding to given electron Lorentz fac- 
tor and to the intensity-averaged synchrotron frequency. Thus 
we expect that the up-scattered radiation is correctly calculated 
at frequencies that are not too close to the limits of the inverse 
Compton spectrum. We checked the correctness of this sup- 
position by calculating burst light-curves using the full shape 
of the inverse Compton spectrum (which leads to substantially 
longer runs) and by comparing them with those obtained using 
the monochromaticity approximations. It is worth stressing here 
that, for an accurate treatment of the inverse Compton losses, 
one has to resort to numerics in order to take into account the 
relativistic beaming of the local synchrotron output, due to the 
relative motion of the cells where the photons are generated and 
up-scattered. 

Light-curves and spectra. We integrate the synchrotron 
and inverse Compton emission over the electron distribution, 
the volume of the entire shocked fluid, and the evolution of 
the interacting shells, to calculate the observed light-curves and 
instantaneous/brightness-averaged spectra, taking into account 
the beaming, Doppler frequency shift and time contraction due 
to the relativistic motion of the radiating fluid. The relativistic 
effects are dependent on the angle between the radial direction 
of outflow and the direction toward the observer. Therefore, 
assuming that the shells are spherically symmetric at least within 
the cone of half-angle T~ 1 visible to the observer, the integral 
over volume is a double one (over shocked cells and over the 
angle relative to the line of sight toward the observer), which 
makes the burst spectrum and light-curve to be a quadruple 
integral. 



4. NUMERICAL SPECTRA AND LIGHT-CURVES 

There are a number of model parameters which affect the 
burst light-curve in the 100 keV-300 keV range spectrum. To 
simplify things, we consider here the interaction between two 
spherically symmetric shells that have fixed Ei = 10 53 ergs, 
£„ = 2 x 10 52 ergs, = 100 and F = 50. The initial kinetic 
energy of the shells can be much lower if they are emitted within 
a relatively narrow cone; these values were chosen so that, for 
a burst located at redshift z = 1 and having the low efficiency 
corresponding to (G = 2, Y = 5), the fluence in the BATSE 
window is above 10~ 7 ergcm~ 2 . To maximize the brightness 
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of the burst we assume that electrons reach equipartition with 
protons and magnetic fields after shock acceleration: e e i = 1/2. 

Equations (A6) and (A7) show that there are only two param- 
eters left that determine the burst spectrum: e mag , on which the 
spectral peaks depend only weakly, and £, the electron injection 
fraction, on which the same peaks have a strong dependence. By 
varying this injection fraction one can study the dependence on it 
of the relative intensity of the synchrotron and inverse Compton 
components, and determine those values of ( that maximize the 
received flux in the BATSE range. Figure 2 shows the shifting of 
the burst emission toward higher energies with decreasing £, due 
to the increase of the electron Lorentz factor. For ( > 10~ 2 the 
inverse Compton emission occurs in the Thomson regime and 
carries most of the burst energy, while for ( < 10~ 3 synchrotron 
dominates over inverse Compton scattering, as the latter takes 
place in the Klein-Nishina regime. The extent in frequency 
of each component is determined by the ratio -jm / 7m of the 
maximum and minimum electron energy of the power-law dis- 
tribution. The7 m isdeterminedbye e i/£andr' s/i i , (eq. [A4]), 
while 7m is set by the details of the electron acceleration and 
its calculation is more ambiguous. An upper limit on it is set by 
requiring that the radiative cooling timescale is longer than the 
shock acceleration timescale, a condition that does not alter the 
shape of the spectrum near its peak unless the injection fraction 
£ is less than 10 -2 . In all other cases shown in Figure 2, we 
have set ^uhm — 100, which gives the correct shape of the 
spectrum at frequencies where most the burst emission lies. 

As shown in Figure 2, there may be more than one component 
that carries a good fraction of the burst emission and each com- 
ponent may extend over a few orders of magnitude in frequency. 
This suggests that the radiation which falls in the BATSE win- 
dow could represent in some cases a rather small fraction of 
the entire emission of the burst. If we also take into account 
that at most half of the available internal energy can be given 
to electrons by shock acceleration, it results that the efficiency 
at which the shells' kinetic energy is transformed into radiation 
visible to BATSE can easily be one order of magnitude lower 
than the efficiency of <~ 10% at which the shocks convert the 
same kinetic energy into internal (for a study of the latter, see 
Kobayashi, Piran & Sari 1997). We can draw the conclusion 
that the overall process that leads to 7-ray emission has an effi- 
ciency of few percent or lower (see also Daigne & Mochkovitch 
1998, Fenimore et al. 1998). Figure 2 also illustrates the fact 
that the radiation detected by BATSE can be either synchrotron 
emission or self-inverse Compton scatterings. The distinction 
between the two cases could be made through the detection 
of simultaneous emission at energies well below or above the 
BATSE range. 

Figure 3 a shows the spectral evolution of a burst whose 
BATSE window emission is due to inverse Compton scatter- 
ings. The large low energy (20 — 50 keV) slope of the spectra 
shown in Figure 3a is a result of the monochromaticity approxi- 
mation described in §3, which was made for numerical reasons. 
Figure 36 shows the spectral evolution of a synchrotron burst. 
The shaded curve represents the average spectrum, calculated 
as an intensity-weighted average of the instantaneous spectra. 
The curve shown with the thick solid line represents a fit with 
the Band function (Band et al. 1993) to the average spectrum, in 
the range 30 keV-3 Me V. The fit is characterized by a = —1.14 
(the low energy index), f3 = —2.38 (the high energy index) and 
E = 326 keV (vF v peaks at (a + 2)E () = 282 keV). 

As shown in Figure 3 c both bursts exhibit a spectral softening 



that can be well approximated as a power-law. It should be 
noted that this is a good fit if the observer would be able to set 
T = when the inner shell is ejected from the burst progenitor. 
Obviously, the observer can time the burst only when it begins, 
in which case the power-law spectral softening would be a bad 
fit. The result shown in Figure 3c should be understood as 
following: there is a T ~ (1/5 - 1/3) T b , with T b being 
the burst duration, such that, if T were measured from To, 
then the spectral softening would be a power-law in T. A 
clearer characterization of the burst softening is illustrated in 
Figure 3d, which shows the spectral peak decays exponentially 
with the photon fluence <t>23 in the middle BATSE channels, 
a feature that was observed in real GRBs by Liang & Kargatis 
( 1 996). Within our model, the softening of the burst is due to two 
factors. One is that the shell fluid shocked later is less dense 
due to the radial expansion, leading to lower internal energy 
densities in the shocked gas, which, for a constant parameter 
£mag, implies lower magnetic fields. Secondly, the radiation 
emitted by the fluid that moves at larger angles off the line of 
sight toward the center of expansion (the "central line of sight") 
is less blue-shifted by the relativistic motion of the source and 
arrives later at observer than the radiation emitted by the shocked 
fluid expanding at smaller angles relative to the central line of 
sight. The Lorentz factor of the shocked fluid does not influence 
the spectral softening, as it is practically constant during the 
two-shell interaction. 

In calculating the spectra shown in Figures 3a and 36, we 
have considered the interaction of two shells in which the ejecta 
are distributed homogeneously. The general spectral softening 
is the same if one considers a shell to be a collection of "mini- 
shells" moving at the same Lorentz factor, i.e. if the two larger 
shells have a layered structure. In this case the burst light- 
curve exhibits pulses associated with each new layer that is 
shocked (sub-pulse structure). The pulses' fluence depend on 
the kinetic energies of the mini-shells, while the separation 
between pulses is determined by the spatial separation between 
layers. The pulse duration, as well as its shape, are determined 
by the angular extension of the region visible to the observer 
(a spherical cap of half-angle opening <~ the thickness of 
the layer and the electron cooling time-scale. 

The light-curve shown in Figure 4a was obtained considering 
an inner group of layers moving at r\ = 100 and an outer one 
moving at r o = 50. The mini-shells in each group have equal 
masses, corresponding to total shell energies Ei — 10 53 ergs 
and E a = 2 x 10 52 ergs, if the ejecta are spherically symmetric. 
The energy release parameters (see figure caption) where chosen 
such that the gamma-ray burst is due to synchrotron emission. 
Figure 4a shows six individual peaks generated by the shocks 
that sweep six outer layers (in this case the inner mini-shells 
radiate mostly outside the BATSE range), and the 50 keV-300 
keV pulse resulting from the addition of these pulses. In Figure 
46 we show the dependence of the shape of the first pulse in the 
burst shown in Figure 4a on the observing energy. The pulse 
lasts longer at lower energy, as observed in real GRBs. If the 

pulse duration is defined by T^ lse = j pulse {F^ / F^ ax )dT, 

where F^ and Fmax are the flux and its maximum value in 
the i-th observing channel, then we find that T^i se 
where Ei is the geometric mean of the upper and lower energy 
limits of channel i. The same dependence is found for the other 
pulses (Figure 4c), as well as in the case when the radiation 
in the BATSE range is due to inverse Compton emission. The 
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exponent changes slightly if the FWHM of the pulse is used 
instead of the integral duration defined above, if one uses the 
lower or upper limits of each BATSE channel instead of their 
geometric means, or if the photon fluxes are used instead of the 
energetic ones. The dependence found by Norris et al. (1996) 
in real GRBs is T„ . oc £!~(°' 3 ~ - 4 ) w hich is stronger than 

pulse i ' t> 

shown by our numerical bursts. 

The usual argument used to explain the observed pulse du- 
ration - energy anti-correlation is that the electron synchrotron 
cooling timescale t sy oc j~ ] oc (hv sy )~ 1 / 2 . Because the syn- 
chrotron spectrum of an electron extends over more than just one 
BATSE channel, the above exponent of 1/2 should be regarded 
only as an upper limit, which is consistent with observations. 
However the argument based on the electron cooling ignores the 
possible contribution of the geometrical curvature of the layer 
and of its thickness to the pulse duration. The spread in photon 
arrival time due to the curvature of the shell is Tg ~ t/(2Y 2 gh ), 
where t is the lab-frame time, because the observer receives 
radiation mainly from the fluid moving within 6 ~ r^ 1 off 
the central line of sight. If the observer frame electron cooling 
timescale T sy ~ t sy /(2Y 2 sh ) exceeds T$, then the shell radius 
r = ct and its volume increase significantly during t sy , lead- 
ing to excessive adiabatic losses and to a lower burst efficiency. 
Thus, an efficient burst is one where T sy < Tg. Our choice of 
model parameters ensures that the electrons are radiative, and 
implies that the pulse duration - energy anti-correlation is not 
due to the electron cooling)! 

In Figure Ad we show the six pulses that form the burst shown 
in Figure 4a, as seen in the 50 keV-300 keV band. The fluxes 
have been normalized to their maxima and the pulses have been 
aligned at their peaks. It can be noticed that these pulses have a 
sharp rise and a slow decay and that are more asymmetric than 
the average pulse shape determined by Norris et al. (1996) for 
separable pulses in long and bright GRBs. The pulse shape is 
determined by the relative importance of the layer thickness, 
its angular spreading and the electron cooling. As mentioned 
before, for our choice of energy release parameters the electron 
cooling is too fast to play any part in determining the pulse 
shape. If the thickness of the layer is not taken into account, 
then the emitting region is approximated by a surface and the 
emission of radiation is almost instantaneous. This explains the 
similarity between the shape of the numerical pulses and that 
obtained analytically by Fenimore, Madras & Nay akshin (1996) 
for the pulse resulting from the interaction between a single 
shell and a stationary medium, when the shell is approximated 
as infinitesimally thin and the radiated power as a delta-function 
in time. The effect of taking into account the shell thickness can 
be assessed by comparing the light-curves shown in Figures 4a 
and 46, corresponding to a shell containing several thin layers 
and one infinitesimally thin layer, respectively. Thus, a shell 
thickness ~ r/Y 2 yields a slightly more symmetric pulse, and a 
shell thickness larger than predicted by Meszaros et al. (1993) 
is required to obtained an even more symmetric pulse. 

The above equation for the spread in the photon arrival time 
due to the geometric curvature of the source (Tg ~ t/(2r 2 sh )) 
implies a correlation between the pulse duration and pulse onset 
time for all the pulses arising from the collision of two inho- 
mogeneous shells. This correlation is illustrated in Figure Ad, 
which clearly shows that later pulses last longer than earlier 
ones. However, such a correlation will not be manifested by all 



the pulses in a real GRB, as it is possible that pulses seen close 
to each other by the observer were emitted by different groups 
of colliding mini-shells, located at different radii, and thus hav- 
ing different durations. This is illustrated in Figure 5, where we 
considered 10 pairs of interacting shells. Each pair has the same 
parameters as the two shells that yield the burst shown in Figure 
4a, except that the time interval t v between their ejection differs 
from pair to pair. Therefore each pair has a different interaction 
radius. The time elapsed between the ejection of successive 
pairs is also considered variable. The light-curve shown in the 
Figure 5a was calculated assuming that, in the frame of the 
shocked fluid, the emission is concentrated within two cones 
of solid angle 47r/5 sr directed outward and inward along the 
radial direction, so that the observer receives radiation from the 
shocked gas moving only within ~ T~ h /2 off the central line 
of sight. This was done in order to reduce the spread in the 
photon arrival time due to the angular spread of the shell and 
the pulse overlapping which is present in Figure 4a. A ran- 
dom injection of several groups of shells can be simulated by 
repeating and superposing the template pulse shown in Figure 
5a, corresponding to arbitrary values of the pulse onset time, 
intensity and duration determined by the value of r, n < for each 
pair of interacting shells. We chose for simplicity a periodic 
shell Lorentz factor and mass distribution in the wind which 
ensures that each shell suffers only one collision. An example 
of such a complex light-curve simulation is shown in Figure 56. 



5. CONCLUSIONS 

As illustrated by the numerical results shown in the previous 
section, the observed GRB 7-rays can be due either to the syn- 
chrotron emission from shocked ejecta, if the electron injection 
fraction is small enough (typically around 10 -2 ) to ensure that 
the accelerated electrons reach high random Lorentz factors, or 
to the up-scattering of the synchrotron photons, if the electron 
injection fraction is not far below unity. The particular choice 
of the shell Lorentz factors used in this work (determined by 
numerical reasons) has lead to an overall efficiency of convert- 
ing the initial shell kinetic energy into 7-rays in the range 25 
keV-1 MeV that is below 1%. A wider range of Lorentz factors 
can increase this efficiency. Nevertheless, since only a fraction 
of the internal energy of the shocked gas can be given to the 
electrons through shock acceleration, and the spectral emission 
range of the burst is very broad, the efficiency in a given in- 
strument's range (such as BATSE) can be substantially smaller 
than that calculated using only the dynamics of the interaction 
between shells. 

The fact that the emission from the gas moving at larger angles 
relative to the central line of sight arrives later to the observer, 
and has a softer spectrum than the emission from the fluid flow- 
ing at smaller angles, leads to an increase in the pulse duration 
with decreasing energy. We find that, if the geometrical curva- 
ture of the shell were the only factor that determines the pulse 
duration, then the pulse duration dependence on energy would 
be E~° , which is weaker than observed. The pulse duration 
is also determined by the electron cooling timescale and by the 
shell thickness. This inconsistency with observations may be 
due to our choice in the numerical calculations of parameters 
that led to an electron cooling timescale smaller than that of 
the adiabatic losses, to maximize the burst efficiency, and to an 



Because the sweeping up of a layer is simulated through the addition of a single cell of shocked fluid (which is done for numerical reasons), it follows that the 
duration dependence with energy shown by the numerically simulated pulses is due only to the geometrical curvature of the emitting shell 
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interaction radius sufficiently large to ensure optical thinness, 
which led to a spread in the photon arrival time due to shell's 
curvature dominating that due to its thickness. 

A softening of the burst spectrum with time is a natural con- 
sequence of the above-mentioned correlation between the angle 
relative to the observer at which the emitting fluid moves, the 
arrival time and the hardness of the radiation received, together 
with the progressive decrease of the turbulent magnetic field in- 
tensity due to the radial expansion of the ejecta. The evolution of 
the break energy of the numerical spectra can be approximated 
quite well as an exponential in the 50 keV-300 keV photon flu- 
ence, for bursts in which either the synchrotron or the inverse 
Compton emission peaks around 100 keV. The spectrum of the 
synchrotron burst is well approximated by the Band function. 

As shown by Daigne & Mochkovitch (1997), a significant 
subset of the spectral-temporal correlations observed can be ex- 
plained within a simple treatment of the kinematics and dynam- 
ics of unsteady winds. Here we have shown that in a more com- 
plete radiation and hydrodynamical treatment, some of the burst 
features found by Daigne & Mochkovitch (1997) are qualita- 
tively confirmed, while some quantitative details differ, possibly 
because of the more detailed physics (including the hydrody- 
namic treatment) introduced here. In particular, our treatment is 
able to reproduce an observational feature that previously was 
not easily obtained, namely a spectral softening in time. 

The "kinematic" light-curves calculated by Kobayashi et 
al. (1997) for internal shocks exhibit a complicated structure, 
and their bolometric temporal profiles bear a good resemblance 
to those of real GRBs. Our band light-curves show a similar 
behavior, and can in addition probe the physical origin of more 
detailed effects, such as spectral-temporal correlations. In the 
case where all pulses within a burst arise from one single group 
of interacting shells, one would expect a correlation between 
the pulse duration and the time measured from the beginning of 
the burst. This is due to the fact that, on average, successive col- 
lisions within a single group of shells take place at larger radii, 
and that the pulse onset time and duration are both proportional 
with the radius where each collision takes place. Obviously, 
the pulses seen in a real GRB may be due to several groups of 



shells interacting at different radii, and thus producing sets of 
pulses of different durations that overlap and mix, as illustrated 
in Figure 5b. 

Whereas a pulse duration increase with time is the signature 
of the ejection of closely bunched shells with different Lorentz 
factors, a lack of continuous correlation between pulse duration 
and pulse onset time would indicate repeated episodes, stretch- 
ing over a longer period of time, of ejection of bunches of shells. 
This may be useful in mapping the injection time-history by the 
central engine, and perhaps shed some light on the dynamics 
of the post-collapse or merger-disruption event. For instance, 
the above correlation (or lack thereof) could be used for testing 
whether double (or multiple) peaked bursts arise from discrete 
and separated "events". Examples of such discrete events might 
be, e.g. , the accretion of discrete rings of disrupted matter (or 
more speculatively, the collapse to a neutron star followed by 
collapse to a black hole or the collapse of a primary followed by 
explosive deleptonization of a small mass neutron star compan- 
ion). Each discrete event would be characterized by the above 
correlation within the event, which then resets itself at the next 
event, if they are truly discrete and independent. 

The results published so far on the GRBs produced by inter- 
nal shocks in an unsteady relativistic wind show that this model 
is able to explain many of the well established properties and 
correlations observed in real bursts. Further work is necessary 
to analyze the model features at a more detailed level: numer- 
ical results having sufficient temporal resolution would allow 
a comparison with the other correlations among pulse features 
(peakedness, asymmetry, width, centroid lag) found by Norris 
et al. (1996). Such features, as well as the general efficiency, 
are issues that may need to be addressed within more specific 
models for the burst progenitor. The advantage of the calcula- 
tions presented here is that they are independent of any specific 
model about the primary event, the only requirement being that 
the central engine produces a sufficiently energetic relativistic 
wind. 

This research was supported in part by NASA NAG5-3801 
and NAG5-2857 
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A. RELATIVISTIC SHOCKS: T/ > Y 

If the shocks are relativistic (G> 1), the solution of equation (Eh leads to 



r' 



sh.i 



G 3 



(G-y)(Gy-l) 



1/2 



sh = 



r,-r, 



' G-y 



1/2 



r' = nV 



(Al) 



provided that 2/G <y < G/2 (otherwise the shock propagating in the denser shell cannot be considered relativistic). From equation 
(Q) and assuming the adiabatic index 7 = 4/3 for a hot gas, the internal energy density in the shocked fluid is 



G 3 y 



[G-y){Gy-l)_ 

and, using equations M) and (0) to derive the pre-shock co-moving densities, the magnetic field can be calculated: 



B - 9 5 x 10V 1/2 R i/2 r -3,-3/2 



G 3 



[G-y){Gy-l) 



1/2 



[G], 



(A2) 



(A3) 



where e mag is a parameter describing the magnetic field strength (see §3) and where the usual notation A = 10 n A n was used. The 
minimum electron Lorentz factor j m of the power-law distribution of shock accelerated electrons is 



(A4) 



where e e i and C parameterize 7,„ and an electron index p = 2.5 was assumed (see §3). Together with equation (Al), equation (A4) 
leads to 



7n 



300 



(G-y)(Gy-l) 



1/2 



7n 



?/ 



The observed peaks of the synchrotron and inverse Compton emission are straightforward to calculate: 



hv sy ^i = 1.6 



J/2 



p l/2 r -2.-3/2 



G 5 



= 140 ( ^ 



p l/2 r _2,-3/2 
fc mag,-l- c 'i,53 i o,2 ( 'i;,0 



4/2 



(G-y) 2 (G2/-l) 
G 8 



[keV] = — 



(G-y) 3 (Gy-l)2 



[MeV] 



hp* 



Y 2 



(A5) 



(A6) 



(A7) 



assuming that the up-scattering of photons takes place in the Thomson regime. 
The lab-frame electron radiative cooling timescale is upper bounded by the synchrotron cooling time 



2.9 



T 



-\ rp-\ r 7 .3 

t mag,-l nj i,53 L o,2 L vfi 



(G-y){Gy-\f 
G A 



]jtsy,o i 



(A8) 



which is much shorter than the timescale for adiabatic losses t ac i ~ ti nt , where, from equation (gj), ti nt = r; m t /c~2x 10 T Q 2 t v fi 
s. Equation (eq. [A8]) shows that for model parameters that are not far from the chosen scaling values, the synchrotron cooling 
timescale of the electrons radiating at hv ~ 100 keV is smaller by four orders of magnitude than the time during which the shock 
sweeps up the shell, which is of order tint. F° r the observer, the electron cooling time appears ~ times shorter, due to the motion 
of the source, and is therefore much smaller than the spread in the photon arrival time due to the shell curvature, which is of order 
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B. SYNCHROTRON AND INVERSE-COMPTON EMISSION 



Synchrotron emission. The calculation of the synchrotron spectrum is based on a numerical approximation derived from the 
equations given by Rybicki & Lightman (1979). The synchrotron power per unit frequency P(u>) (for one electron) is 



3 5/2 p 



UJ 



P( U ) = -—--2LF[ — ); (Bl) 

where P sy = (1 /6n)aThcB 2 (^l — 1) is the frequency-integrated synchrotron power (and averaged over the pitch angle), uth 
being the cross-section for electron scattering and j e the electron Lorentz factor, and lu c = (3ir /8) (eB / 'm e c)^j 2 is the synchrotron 
frequency (averaged over the pitch angle), with e and m e the electron charge and mass, respectively. F(u) can be approximated by 
F(u) ~ 1.78 u 0,297 e~" for 10~ 35 < u < 10 0,5 with a maximum error of 5%. For u such that F(u) > 0.5 (i.e. close to the peak of 
the synchrotron spectrum) this approximation is accurate to better than 1%. At frequencies far from the synchrotron peak, we used 
the approximations given by Rybicki & Lightman (1979): F(u) ~ 2.15 u 1/3 for u <C 1 and F(u) ~ 1.25 e~ u ^/u for u > 1. 

Inverse Compton scatterings. The peak of the inverse Compton spectrum and the inverse Compton losses are calculated using 
approximations derived from the equations given by Blumenthal & Gould (1970). The inverse Compton spectrum peaks at the energy 
e p given by 

_ 4g p 7 2 e 



1 + (^q P le£o/m e c 2 ) 



(B2) 



where eo is the energy of the incident photon, and q p is a factor that depends weakly on j e eo. We found that q p can be approximated 
with an error below 1 % by 

_ 1 5.91 x 10~ 2 5.09 x 10~ 2 
qp ~ 2 + 1 +0.184W 1 - 31 + 1 +51.6U 1 - 45 ; (B3) 

with v = 7 e eo/m e c 2 . Equations (B2) and (B3) lead to q p = 0.610 and e p ~ 2.447 2 eo in the Thomson regime (v <C 1), and to 
q p = 1/2 and e p ~ j e m e c 2 in the extreme Klein-Nishina regime (v ^> 1). For the inverse Compton power P ic (per electron) we 
found that the following approximation: 

Pic(v) / [1+7.67 exp(2.431ogi;)]- 1 v<\ 

P^h ~ \ 0.107w- 107 exp(-0.5691og 2 i;) 1< v < 10 L5 ' K ' 

where P? h = (4-/3)<7ThcU 2 y (jl — 1) is the inverse Compton power in the Thomson regime, and U sy the energy density of the 
photon field that is up-scattered, has a relative error that increases with v, reaching a maximum value of 10% at v = 10 L5 , where 
Pic ~ 10~ 3 P? h . The inverse Compton losses are severely reduced by the Klein-Nishina effect at v > 10 1 5 and an accurate treatment 
of these losses is not necessary. 
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Fig. 1 . — Left panel: Lorentz factors of the shocked fluid as measured in the frames of the yet un-shocked gas in the inner shell (F' sh i 
- thin lines) and in the outer one ( F' sh - thick lines), for a range of shells' density ratio and for three values of the ratio of their 
lab-frame Lorentz factors. For given ratio rj/r o , the Lorentz factors F' sh • and F' sh are symmetric relative to the ordinate. Right 
panel: the efficiency of the shocks in converting the shells' total kinetic energy into internal energy. The legend is the same as for the 
left panel. 



10 



w 10° 



E 
o 

CD 



10° 



10" 



to 10 

E 

o 

o5 10 



10"' 



; 1 SY rs 

S~ \ IC rs-rs - 

t- / ^ SYfs -= 

; i ^ ic fs-fs - 

A\ 1 ^ 

r [\ | 1 | | j 


: C=0-1 

; i i 




: £=0.01 

; / A M 

/ J\ : 


: £=0.001 

/ \ 

r /-■/-.. \ ] 
: / / \ : 

/ / W /"■. - 

z / i \\ \ \ 

- / / \1 / V 



10 10 10 10 

hv [eV] 



10 1 



10 J 



10 5 10 7 

hv [eV] 



10 a 



= 10* 



10" b 



10" 



CD 

— \ 

CQ 
O 

3, 

if) 



-- 10" c 



10"" 



4 10" 



CD 
CQ 

o 

3, 

rv 



10 1 



Fig. 2. — Dependence of burst spectra on (, the electron injection fraction, for fixed Ei — 10 53 ergs and E Q = 2 x 10 52 ergs in 
4n sr, Fi = 100, r o = 50, and t v = 1 s, e e j = 1/2, e mag = 10 _1 and p = 2.5. The source is located at redshift z = 1, with 
Hq = 75kms _1 Mpc _1 and £2 = 1. The synchrotron (thin lines) and self-inverse Compton (thick lines) emissions from the reverse 
shock (rs) and forward shock (fs) are shown separately. The mixed components resulting from up-scattering by one of the shocks of 
synchrotron photons generated by the other shock are too weak and do not appear in the graphs. Note that for ( <; 10 _1 , the observer 
receives up-scattered photons in the BATSE range, while for ( £ 10~ 2 the 7-ray burst is due to synchrotron emission, and the inverse 
Compton emission is diminished by the Klein-Nishina reduction. 
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Fig. 3. — Spectral evolution of bursts with same parameters as in Figure 2, except £ = 0.5 for panel a (inverse Compton emission) and 
p = 3 for the panel b (synchrotron emission). The average spectrum shown in graph b with a shaded continuous line is fit in the range 
30 keV-3 MeV by the Band function with parameters a = — 1.14, /3 = —2.38 and Eq = 326 keV. Graph c shows in log-log scale 
the decrease with time of E p , the peak of vF v (power-per-decade), i.e. the softening of the burst spectrum, x's are for synchrotron 
emission and pluses for inverse Compton radiation in BATSE window. For both types of bursts, a power-law approximates quite well 
the spectral softening. Graph d shows the exponential decay of the peak E p with the photon flux 4>23 in the range 50 keV-300 keV. 
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Fig. 4.— 

a : Light-curve of a burst arising from the collision between an inner set of faster and more massive shells with a group of six outer, 
slower and less massive shells. The Lorentz factor and total kinetic energy of the shells is each group are = 100, Ei = 10 53 ergs 
(assuming spherical symmetry), and r o = 50, E a = 2 x 10 52 ergs. The burst is located at redshift z = 1. Other parameters are 

s el = 1/2, C = lO" 2 , s mag = lO" 1 . 

6 : The first pulse in graph a, as seen in each BATSE channel. The pulse lasts longer at lower energies, 
c : The dependence on observing energy of the duration of the first three pulses shown in graph a. 

d : Evolution of the pulse duration with pulse onset time. For a collision between two layered shells, later pulses last longer than 
earlier ones. 
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Fig. 5. — Light-curves arising from pairs of shells having the same parameters as the pair that yields the burst shown in Figure 4a, for 
a pre-beaming factor 4-ir /5 of the comoving radiation. Panel a shows the pulse from a single pair, while panel b shows the light-curve 
from a set of 10 pairs (see text for details). 



